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Abstract. EK Dra (HD 129333) is a young, active, nearby star that is orbited by a low mass companion. By 
combining new speckle observations with old and new radial velocity measurements we find that the orbit is highly 
eccentric with e = 0.82 ± 0.03, and we derive the true masses of both components. The masses are 0.9 ± 0.1 M© 
and 0.5 ± 0.1 Mq, for the primary and secondary, respectively. From high resolution spectra we derive a new Tcfi 
of 5700 ± 70 K and a log g of 4.37 ± 0.10, which is different to previous estimates. However, the new spectroscopic 
distance differs by only 5.8% to the distance derived by parallax measurement by the Hipparcos satellite and thus 
the stellar parameters are presumably more realistic than older determinations. We derive a somewhat higher 
value for the metallicity of [Fe/H] = —0.16 ±0.07. EK Dra turns out to be one of the few nearby young stars that 
will evolve similarly to the Sun. The precise radial velocity measurements taken in the course of this program also 
allow us to shed more light on the activity of this star. In 2001 and 2002 we find radial velocity variation with a 
period of 2.767 ± 0.005 days which we interpret as the rotation period. This signal vanishes in 2003. However the 
signal can be recovered if only the spectra in which the photospheric lines are asymmetric are used. On the other 
hand, we do not find a close correlation between the asymmetry of photospheric lines and the radial velocity. 

Key words. Stars: individual: EK Dra, activity, fundamental parameters, binary: spectroscopic and visual 



1. Introduction 

EK Dra (GJ 559. lA, HD 129333) is a star that has roughly 
the mass of the Sun. The Henry Draper Catalogue and 
Extension (Cannon & Pickering 1918-1994) and Simbad 
■ database list a spectral type of F8 for this star. The equiv- 
' alent width of Li I is about 0.2 A, and therefore it must 
. be relatively young. Frohlich et al. (2002) and references 
' therein discuss the age and activity connection. 

Various studies of EK Dra at different wavelength 
regimes have been carried out including longterm photo- 
metric monitoring over decades. Dorren & Guinan (1994) 
have observed strong variable chromospheric emission 
lines in their UV spectra. The star has the highest known 
CaH H and K emission level of any known early G-type 
star which is not a close binary (Soderblom 1985). 

The star is rapidly rotating (w sini — 16.5 ± l.Okm/s) 
and has dominant spot features at ^ 70° — 80° that 
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could also be the offshoot of a large polar spot. These 
spots are located at a higher latitude than typical spots 
on the Sun (Strassmeier & Rice 1998). Strassmeier & 
Rice (1998) measure several rotation periods between 
2.599 ± 0.001 days and 2.796 ± 0.026 days using different 
methods where for their purpose they adopt a longterm 
photometric period of 2.605 days. Coronal emission was 
also observed in X-rays and as well in the radio regime. 
The X-ray light curve is significantly variable, with the 
emission from the cooler plasma being strongly modulated 
by the rotation period, while the emission from the hotter 
plasma is only weakly variable (Guedel et al. 1995). 

A 12 to 14 year cyclic variability was discovered by 
Dorren & Guinan (1994) and Dorren et al. (1995) using 
photometric and spectroscopic data. They observed that 
the Ca n H & K emission index increased during that time. 
A decline of brightness since 1994 was noted by Frohlich 
el al. (2002). The star became fainter as its mean level of 
chromospheric activity rose. These findings can be inter- 
preted as signs of a spot cycle. 

EK Dra also is a long period binary star where the 
secondary is much fainter than the primary. Duquennoy 
& Mayor (1991) used a period of 11.5 years in order to 
derive the first preliminary spectroscopic orbit. However, 
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as will be discussed in the next section, the true period is 
45 ± 5 years. 

EK Dra is a wcU-studicd young, active and nearby star. 
It thus serves as one of the best-studied young stars evolv- 
ing similar to the Sun. Since it is a long period binary, it is 
possibly one of the few cases of young stars for which the 
true masses can be determined. The aim of this paper is to 
derive the true mass, and to calculate a new atmospheric 
model. These will then allow us to compare the proper- 
ties of this object with evolutionary tracks. Additionally, 
we have obtained a large number of radial velocity mea- 
surements which will give us new insights into the stellar 
activity and the influence of stellar activity on precise ra- 
dial velocity measurements. 

2. Deriving the true masses of EK Dra A and 
EK Dra B 

By combining the data from our speckle interferometry 
and RV data from the literature with our own RV mea- 
surements, it is for the first time possible to derive the 
true masses of EKDraA and EKDraB. 

2.1. Speckle interferometry : observations 

EK Dra was repeatedly observed with the 3.5m-telescope 

on Calar Alto from 1991 to 2002 using speckle interfer- 
ometry, mostly in the K band (A = 2.2 /xm). An overview 
of these observations is given in Table 1 and is displayed 
in Fig. 1. Observations 1 and 2 made use of a device for 
one-dimensional speckle interferometry. For these obser- 
vations, the technique of data acquisition and reduction 
has been described by Lcincrt & Haas (1989). All other 
data points were obtained using the near-infrared cameras 
MAGIC and OMEGA Cass that are capable of taking se- 
quences of short exposures (typically r « 0.1 s), and thus 
allow speckle interferometry with two-dimensional detec- 
tor arrays. Details of data reduction and analysis have 
been described by K5hler et al. (2000). Briefly, 1000-1500 
short exposures ( "frames" ) are taken for EK Dra and the 
nearby point source (PSF calibrator) BS 5436. BS 5436 
is an F4IV star at a distance of 31.5 pc which implies 
that its relative size is < 0.3 mas, and thus it is consid- 
erably smaller than the resolution of 1.22A/_D = 160 mas 
of the 3.5-m-telescope in the K-band. In Fourier space 
this reduces to \/D — 130 mas. These images are stored 
in "data cubes" of 250 frames. The telescope position is 
switched between object and PSF calibrator after each 
data cube to observe both under nearly identical atmo- 
spheric conditions. After background subtraction, flat- 
fielding and correcting for bad pixels, the data cubes are 
Fourier transformed. The modulus of the complex visi- 
bility is derived by deconvolving the power spectrum of 
the object with that of the PSF calibrator, while the 
phase is reconstructed using the Knox-Thompson algo- 
rithm (Knox & Thompson 1974) and also the bispectrum 
method (Lohmann et al. 1983). The complex visibility is 
averaged over all observations of EK Dra taken on one 
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Fig. 1. Orbital motion of the companion EK Dra B from 
1991 to 2002. The primary is located at (0,0). As can be de- 
duced from Table 1., the last two points are separated by one 
year. Primary and secondary are thus close to the apastron. 
Combined with the RV measurements we estimate a period of 
about 45 ± 5 years because Duquennoy & Mayor observed the 
RV minimum in May 1987. 



night. Finally, the binary parameters - position angle, 
projected separation and flux ratio Fb/Fa - are derived 
from a model flt to the complex visibility in Fourier space. 
Except for the data points 1 and 2 in Table 1, for all ob- 
servations the relative astrometry of the components has 
been put into a consistent reference frame. This reference 
frame is primarily based on astrometric flts to images of 
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Table 1. An overview of all spatially resolved observations of the binary system EK Dra AB, carried out with near-infrared 
speckle interferometry at the 3.5-m telescope on Calar Alto. 



No. Date 


Epoch 


Position 


Projected 


Filter Flux ratio 


Instrument 






angle [°] 


separation [mas] 


Fb/Fa 





1 


19.03.1991 


1991.2135 


173.9 ± 4.3 


282 ± 12 


K 




ID 


2 


13.02.1992 


1992.1232 


171.7 ± 3.9 


313 ± 13 


K 




ID 


3 


05.10.1993 


1993.7611 


173.7 ± 0.9 


456 ± 14 


K 


0.066 ± 0.003 


MAGIC 


4 


26.01.1994 


1994.0712 


174.0 ± 0.3 


486 ± 10 


K 


0.084 ± 0.003 


MAGIC 


5 


13.12.1994 


1994.9501 


175.5 ± 0.4 


501 ± 5 


H 


0.087 ± 0.003 


MAGIC 


6 


22.11.1997 


1997.8926 


174.8 ± 0.6 


647 ± 12 


K 


0.084 ± 0.004 


MAGIC 


7 


10.02.2001 


2001.1123 


172.8 ± 0.6 


679 ± 5 


K 


0.111 ± 0.006 


OMEGA Cass 


8 


03.11.2001 


2001.8406 


174.0 ± 1.0 


709 ± 10 


K 


0.081 ± 0.008 


OMEGA Cass 


9 


21.10.2002 


2002.8049 


171.9 ± 1.0 


725 ± 8 


K 


0.081 ± 0.010 


OMEGA Cass 



the Orion Trapezium cluster core, where precise astrome- 
try has been given by McCaughrean & Stauffer (1994). 

2.2. Radial velocity measurements: data 



• this paper 

o Duquennoj & Mayor 1991 



HJD-345OO0O 



Fig. 2. The RV data from Duquennoy & Mayor (1991) marked 
with open circles and the three years of RV measurements ob- 
tained at the TLS marked with filled squares. 



at Kitt Peak. The latter enables us to compute the relative 
velocity shift between stellar and iodine absorption lines as 
well as to model the temporal and spatial variations of the 
instrumental profile; see Valenti et al. (1995) and Butler 
et al. (1996) for a description of the principles behind this 
technique. Fig. 2 shows our RV measurements together 
with those obtained by Duquennoy & Mayor (1991). 

RV measurements have been made at TLS since 2001 
and these show that we can achieve a routine RV preci- 
sion of « 3ms~^. However, our RV measurements for EK 
Dra have an error of about 30ms~^. Two factors degrade 
the RV precision of our EK Dra measurements. First, EK 
Dra has a vsmi of 16.5 ± l.Okm/s. Since the RV error is 
proportional to the vsmi, the error compared to a more 
slowly rotating star with comparable S/N should be sev- 
eral times worse. Second, EK Dra is an active star and as 
demonstrated by Saar & Donahue (1997) the activity can 
introduce significant RV "jitter" depending on the level of 
activity of up to several tensof m/s. 

The relative RV measurements obtained by us were 
converted to absolute values by measuring the absolute 
RV of the template spectrum by fitting Voigt functions to 
photosphcric lines with an equivalent width larger than 
0.1 A in the wavelength range from 5000 to 6000 A. 



EK Dra is one of the stars monitored during the RV 
search program for young and active stars of the Thiiringer 
Landessternwarte (TLS) described by Hatzes et al. (2003). 
For this program we use the 2-m- Alfred Jensch telescope 
of the TLS, which is equipped with an echelle spectro- 
graph with a resolving power of AA/A = 67000. During 
the observations an iodine absorption cell is placed in the 
optical light path in front of the spectrograph slit. The 
resulting iodine absorption spectrum is then superposed 
on the stellar spectrum providing a stable wavelength ref- 
erence against which the stellar RV are measured. In the 
first step, the spectra are bias-subtracted, flat-fielded and 
extracted using standard IRAF routines. 

In the second step the RVs are calculated by model- 
ing the observed spectra with a high signal-to-noise ratio 
template of the star (without iodine) and a scan of our 
iodine cell taken at very high resolution with the Fourier 
Transform Spectrometer of the McMath-Pierce telescope 



2.3. Deriving the masses of the two components 

The motion of EK Dra B with respect to the primary 
from 1991 to 2002 is shown in Fig. 1. Although almost 
no curvature is seen, this relative motion is definitely not 
caused by the proper motion of EK Dra A with respect to 
a chance-projected background star. The proper motion of 
EK Dra A is = — 138.6mas/yr and = — 11.9mas/yr 
{Hipparcos catalog), which is roughly perpendicular to the 
observed motion and also much faster. Furthermore, the 
companion has significantly slowed down over the time- 
span covered by the observations, which can only be ex- 
plained by orbital motion. The slow-down additionally im- 
plies that the last observations are close to the apastron. 
The spectroscopic observations of Duquennoy & Mayor 
(1991) cover the periastron. By combining the radial ve- 
locity data with the speckle observation it is possible to 
derive a first orbit of the system and the true masses 
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Fig. 3. The projected distance versus time. The projected dis- 
tance was converted in to AU by using the Hipparcos distance. 
The line shows the fit using the elements given in Table 2. 




Fig. 4. The radial-velocity curve together with the fit using 
the values given in in Table 2. 

of both components using the distance as measured by 
Hipparcos. The projected distance versus time is shown 
in Fig. 3. The deceleration is clearly visible. By combin- 
ing the speckle imaging data with the radial- velocity data 
(Fig. 4), we derive a period of 45 ± 5 years. Given the or- 
bital period the radial velocity (Fig. 4) constrains very well 
the mass-ratio of the two components. By fitting an orbit 
to th{^ speckle data (Fig. 5) combined with the information 
on the projected velocity and the Hipparcos distance, we 
can calculate the true masses of the components and all 
other orbital elements by solving the Keplerian equations 
(Kepler 1609, Kepler 1618). The orbital elements are sum- 
marized in Table 2. For the masses of the two components, 
we find 0.9 ±0.1 Mq and 0.5 ±0.1 Mq, for the primary and 
the secondary, respectively. With e = 0.82±0.03, the orbit 
is surprisingly eccentric. The distance between the stars 
is nevertheless still 2.2 AU at the periastron. It thus seem 
unlikely that the secondary has a big impact on the activ- 
ity level of the primary. 



3. Spectral synthesis analysis 

We observed EK Dra on September 11, 2001 and May 
25, 2002 from Calar Alto using the high resolution echelle 



Fig. 5. Using the distance as measured by Hipparcos and the 
radial velocity measurements it is possible to derive the true 
masses of both components. Shown here is the derived orbit 
using the values given in Table 2. 

Table 2. Orbital elements 



element 



value 



P 

To [HJD] 
e 

Vo 
i 

dist 

f{m) 

Ml 

M2 



(mi-Hmg)^ 



45 ± 5 yrs 
246960 ± 30 
0.82 ± 0.03 
-21.2 ±0.2 km/s 
5.8 ± 0.2 km/s 
11.2 km/s 
14 ± 0.5 AU 
85 ±3° 
-98 ± 1° 
180 ± 10° 
33.94 ± 0.72 pc 

0.06 ± 0.02 M0 

0.9 ±0.1 M0 
0.5 ±0.1 Me 



esitimated value using the mass-ratio of the primary and 

secondary, and Ki 

semi-major-axis of the true orbit 

position angle of the line of intersection between the tan- 
gential plane of projection and the projected orbital plane 
periastron longitude 



spectrograph FOCES (Pfeiffer et al. 1998) mounted on the 
2.2 m telescope. Data reduction and analysis were carried 
out using the reduction pipeline written in IDL especially 
for this fiber-coupled spectrograph. 

For the spectral synthesis analysis we used the model 
atmosphere code MAFAGS. For a detailed description 
of the methods see Fuhrmann et al. (1997). As de- 
scribed there, we deduce the effective temperature from 
the Balmer line wings and the surface gravity from the iron 
ionisation equilibrium and the wings of the Mg lb lines. 
The analysis is performed strictly relative to the Sun. The 
method for determining all stellar parameters was tested 
and compared extensively in Fuhrmann (2004). 
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name T^tt 


logs 


[Fe/H] 


vsini Mv Mboi 


EK Dra 5700 ± 70 


4.37 ±0.10 


-0.16 ±0.07 


16.50 ± 1.00 7.60 4.79 
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Fig. 6. The lithium absorption line at 6707.8 A in the spectrum of EK Dra compared to reflected sunlight on the moon. Note 
the high vsini of 16.5 ± l.Okm/s of EK Dra. The strong lithium line indicates that EK Dra is young. 



3.1. Results for EK Dra 

The chromospheric activity, the variability, the presence 
of core filling-in of the Ha, the calcium H & K and the 
magnesium Ib-lines and a strong lithium absorption line 
at 6707 A (Fig. 6) indicate that the star is indeed young 
and if we beheve the star to belong to the Pleiades, we 
can assume an age of about 125 Myr. 

The spectral syntheses analysis we performed was chal- 
lenging because it is known that the star rotates fast 
(16.5 ± Ikm/s) and it exhibits huge stellar spots which 
are much cooler that the surrounding surface. Depending 
on the position of the star spot, the star could appear 
cooler than it actually is. On the other hand the Ha-line 
core is filled-in up to a level of 0.5. The H/3-line is not 
noticeable fillcd-in. To measure the effective temperature 
we calculate a grid of line profiles of the Ha- and H/3-lines 
using the surface gravity, and iron abundance determined 
using the Fe I & H- and Mg Ib-lines. We fit the wings 
of the strong lines but not the core. The final effective 
temperature is measured using the Ha-temperature and 
comparing it with the H/?-temperature giving a weight of 
75% to the H/?-temperature. If the new derived tempera- 



ture changed more than 50 K compared to the previously 
obtained temperature we repeat the determination of the 
iron abundance as well as the surface gravity. 

In the case of the spectrum taken on Sep. 11, 2001, 
both temperatures were identical while in the spectrum 
taken on May 25, 2002 the H/3-line leads to a 40 K hotter 
star. For a fit to the Ha-line profile see Fig. 7. However, 
the temperature most consistent with all measurements is 
5700 ± 70 K. We have double-checked the derived effective 
temperature with a fit to the Ha and H/3 lines, assuming a 
spot temperature of 4500 K and a spot coverage of 1/4 of 
the visible surface. The resulting measured temperature 
does not significantly change the previously obtained re- 
sults of the one-temperature fit because the continuum of 
the spot is only 10% of the continuum level of the sur- 
rounding stellar surface. This is because the spot is much 
cooler than the photosphere, and hence its total contribu- 
tion to the light emitted from the star is very small. This 
results in a correspondingly small contribution to the total 
spectrum. In fact, the spot only produces a little hump in 
the photospheric line-profiles used for the Dopplcr imag- 
ing, which changes the EW of a photospheric line only by 
a small amount, and cannot be seen in Ha and H/3. Thus, 
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Fig. 7. The Ha-hne together with two theoretical Ha-line profiles to estimate the temperature. The real temperature is in 
between the plotted ones. Note the filling in of the Ha-line. 



even for a highly spotted star, the average line spectrum 
is dominated by the photosphere. 

Comparing the spectroscopic distance determined by 
us to the Hipparcos parallax, wc have a 5.8% discrepancy. 
This makes us confident that the spectral parameters we 
derived, especially the surface gravity logg, are reason- 
able. Regarding the Ha line depth and the line wings 
(Fig. 11, lower panel), we see significant changes of the 
effective temperature from 5700 K to 5580 K when using 
only the Ha temperature. The H/3 temperature is not that 
strongly affected. 

The analysis from Wyse & Gilmore (1995) lead to a 
soniciwhat lower iron abundance of —0.214 which was es- 
timated by narrow-band Stromgen photometry. The tem- 
perature of EK Dra of 5930 K estimated by Dorren & 
Guinan (1994) is higher than our measured temperature of 
5700 ± 70 K. But in our case the Ha- and H/3-line profiles 
would not support such a high temperature in our spec- 
trum. Eggen (1998) also estimated a metallicity of EK Dra 
of — 0.24dex using Stromgen photometry, which is lower 
than the metallicity measured by us of — 0.16dex. 

For their Doppler imaging, Strassmeier & Rice (1998) 
need stellar parameters as an input and they used Kurucz 
model atmospheres where they fixed the iron abundance 
to [Fe/H] = 0.0 and the surface gravity to \ogg = 4.5 to 
derive T^s = 5870 ± 50 K and sin i = 17.3 km/s. They re- 
cover stellar spots at high latitude with a temperature dif- 
ference of AT = 400 K. A polar spot cannot be confirmed 
or excluded by their data. The inclination of the stellar 
rotation axis is ~ 60°. The inclination of the orbit thus is 



significantly different from ^ 60°. Because of the signif- 
icant difference from our spectral parameters (especially 
the effective temperature and the iron abundance) we pro- 
pose; caution when using the conclusions of the Strassmeier 
& Rice (1998) Doppler imaging. 



4. Space motion and age of EK Dra 

Soderblom & Clements (1987) 

estimate the age of EK Dra to be ~ 70 Myr using activ- 
ity indicators and claim that the star could have trav- 
eled from the Pleiades to the solar vicinity with a pe- 
culiar velocity of only 2 km/s. On the other hand, based 
on activity indicators Wichmann et al. (2003) claim that 
the star is even younger than the Pleiades with an age 
of ~ 50 Myr and call it a member of the local associa- 
tion. Wichmann & Schmitt (2003) have traced the space 
motion backward in time and exclude a former member- 
ship of the young associations Lupus-Centaurus-Crux or 
Upper-Centaurus-Lupus. 

Stauffer et al. (1998) derived an age of the Pleiades 
of about 125 Myr using the aproproate distance scale 
and modern stellar evolution calculations. Estimates by 
Basri, Marcy & Graham (1996) using two brown dwarf 
members of the Pleiades give an age of 115 to 125 Myr. 

We have used the proper motion /Xq = —138.61 ± 
0.72 mas, fig — —11.92 ± 0.56 mas and the parallax tt = 
29.46 ± 0.61 measured by Hipparcos, as well as the RV of 
—23.1 ±0.2 km/s measured by us to calculate the galactic 
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space motion {U/V/W) = (1.3 ± 0.3/24.8 ± 0.4/ - 0.3 ± 
0.3)kni/s). 

EK Dra is located in the vicinity of the Sun at a dis- 
tance measured by Hipparcos of 33.94 ± 0.72 pc. EK Dra 
is a young object and it is likely that its age lies within 50 
to 125Myrs depending on the criteria one applies, e.g. if 
one only regards the activity (50Myr) or if one assumes 
it is a Pleiades field star (125 Myr). 

With the average absolute brightness of the primary 
oiMx = 3.4 ±0.1 mag, and ofM^ = 6.0 ±0.2 mag for the 
secondary component, and using the evolutionary tracks 
pubhshed by Baraffe ct al. (1998) (Model: [M/H]=0, a = 
1.9, Y=0.282), one derives a mass of 1.025 ± 0.100 M© 
and an age of 35 Myr (lower limit: 30 Myr and upper limit 
1.5 Gyr on the main sequence (MS)) for the primary. The 
mass derived from the evolutionary tracks agrees with the 
true mass within the errors. The companion is expected 
to be about 6 mag fainter than the primary in the V-band. 
The visible spectrum of EK Dra is thus completely domi- 
nated by the primary. 

5. Activity, line asymmetry and RV variations 
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Fig. 9. The RV data phase folded with the period of 2.769 days 
and split between the different years of observation. The filled 
symbols represent the asymmetric lines and the open circles 
the symmetric lines. 
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Fig. 8. Top panel the bisector analysis of the chromospheric 
inactive lines outside the region which is affected by the iodine 
Hues. Lower panel the original RV data observed in the years 
2001, 2002 and 2003 at the TLS. 

Because EKDraB was close to the apastron at the 
time of the RV observations, and since the period of this 
highly eccentric orbit is 45 ± 5 years, the influence of the 
secondary on this section of the RV curve is negligible 
(Fig. 4). With an accuracy of our relative RV measure- 
ments of 30 m/s it is possible to study the RV variations 
caused by stellar activity. Saar & Donahue (1997) have 
studied the activity-induced RV variations and found a 
relation of As « 6.5/5^1; sin j, where fs is the star-spot 
area coverage [fs is in %, vsmi in km/s, and As in m/s]. 



Fig. 10. The results of the bisector analysis for all three years 
of observations. The filled squares are the data from 2001, the 
crosses the data from 2002, and the circles from 2003. The 
spectra from 2001 mostly have asymmetric lines, whereas in 
2002 we observe asymmetric as well as symmetric lines. In 2003 
the lines are mostly symmetric. The dashed line indicates the 
border between the symmetric and asymmetric spectral lines. 



For spots located at the equator, the filling factor is given 
by this equation. For spots located close to the pole, the 
RV variations would give a rise to a too small filling factor. 
Because spots on the pole do not give rise to RV varia- 
tions, the amplitude of the RV variations in principle just 
give a lower limit for the filling factor fs- From the mea- 
sured RV variations we thus derive (Fig. 8, upper panel) 
the peak-to-peak amplitude variation of about 550 m/s. 
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Fig. 11. The upper panel shows the average profile of Ha of all 
spectra where the photospheric lines are symmetric in compar- 
ison to the average profile of all spectra where the photospheric 
lines are asyrmnctric. The difference between the two Ha pro- 
files is remarkably small. The only notable difference is that the 
core of Ha is a little deeper in the case of the spectra where 
the photospheric lines arc asymmetric. The lower panel shows 
that there are significant variations in the depth of Ha. Also 
the wings of the Ha lines show a difference and lead to a 120 K 
cooler star. This effect is smaller in H/3. 

which impUes that fs > 5.2% clue to the fact that we 
suspect the presence of a polar spot. 

Performing a period analysis wc find a period of the 
RV variation of 2.767 ± 0.005 days in the first two years of 
observations but adding the third year the period search 
programs do not recover the period clearly. When we in- 
spect the data (sec Fig. 9), we see that the amplitude of 
the RV signal declines. 

Moreover, wc sec a change in the bisector of the spec- 
tral lines over the three years of observations. We use the 
classical Doppler imaging hues, Fel at 6400 A, and Cal 
at 6439 A to measure the amplitude of the asymmetry be- 
cause they arc not blended by other weaker spectral lines 
and there are no iodine lines in this part of the spectrum. 
Shown in Fig. 8 (upper panel) is the average velocity dif- 
ference between the lower and the upper half of the spec- 
tral lines. A negative sign means that the lower part is 
blue-shifted with respect to the upper part of the line. 
During the three years the behavior of the bisector sig- 
nificantly changed. In 2001, the velocity difference is rel- 
atively constant with a shift of 700 to 1500 m/s, whereas 
in 2003, the bisector amplitude varies between 1200 m/s 
and —100 m/s. We then can split the lines between asym- 
metric lines with a bisector from 700 m/s to 1500 m/s and 
symmetric lines with a bisector of 600 m/s to —500 m/s 
(Fig. 10). Excluding the symmetric lines from the pe- 
riod analysis we recover the period over all three years. 
Compared to the amplitude of the RV signal of the sym- 



metric lines, the signal for the asymmetric lines is more 
pronounced (Fig. 9). 

Wc can summarize our finding as follows: in 2001 and 
2002, we observe a periodic RV signal and the spectral 
lines were quite asymmetric. This RV signal is less obvious 
in the data from 2003 but can be recovered if only the 
asymmetric lines are used. Clearly, this is the signal of 
stellar spots not located at the pole; the Doppler imaging 
by Strassmeier & Rice (1998) shows the presence of spots 
close to the equator. One possible explanation would be 
that in 2003, there were fewer spots at mid latitudes, or 
alternatively, there were so many plage regions so that it is 
difficult to recover the periodic RV signal from the spots. 
Plage regions are known to produce less asymmetric lines. 

In order to find out if changes in the symmetry of the 
lines are related to the activity level, we averaged all spec- 
tra with symmetric photospheric lines and all spectra with 
asymmetric photospheric lines. Fig. 11 shows the Ha- line 
profiles of all spectra where the photospheric lines are sym- 
metric in comparison with the average profile of all spec- 
tra where the photospheric lines are asymmetric. The dif- 
ference between the two Ha profiles is remarkably small. 
However, we do see significant variations of the depth of 
Ha. Thus testing if the amplitude of the RV signal is de- 
pendent on the depth of Ha, we calculate the average RV 
of all spectra where Ha is deepest and all spectra where 
Ha is shallowest. We do not sec any significant difference 
between the two samples, as the average RV of the first 
sample is —40 ± 80 m/s°, and the average of the second 
— 130 ± 120 m/s". There is also no significant difference 
in the asymmetry of the two, as the average asymmetry 
of the first sample is 860 ± 380 m/s° and the asymmetry 
of the second sample is 1000 ± 350 m/s°. Thus, there is 
no obvious correlation between the RV or the asymmetry 
and the depth of the Ha-line. Additionally, there is only 
a very weak correlation between the asymmetry and the 
RV. 

The change of the temperature as seen in Fig. 11 (lower 

panel) of about 120 K measured by the Ha-line profile is 
not as visible in the H/?-line profile. There, it only shows 
a temperature change of 20 K which lies within the errors 
of this method. Giving more weight (75%) to the H/3 tem- 
perature will result only in a small temperature change. 
But of course the temperature change will also affect the 
Fc I and Fc II ionisation equilibrium and this will have 
an effect on the determination of the surface gravity logg. 
We cannot verify this here because the spectra obtained at 
TLS are mostly observed with the iodine cell in the light 
path which superposes iodine absorption lines on the part 
of the spectrum containing the necessary iron lines. 

How can wc interpret these findings? A spot that is lo- 
cated close to the equator would lead by its appearance to 
a positive asymmetry and also to a redshift of the spectral 
line. When the spot is receding due to rotation, it would 
lead to a negative asymmetry and to a blue-shift of the 
line. In this case we would expect a clear correlation be- 
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tween the asymmetry and the RV. Such a correlation is 
for example observed in RS CVn systems (Donati et al. 
1994). This is clearly not what wc observe. The Dopplcr 
imaging of this star mainly shows a polar spot, and our 
estimate of the spot size from the amplitude of the RV 
variations indicates that there are only a few spots close 
to the equator. 

In the case of sunspots it is well known that the 
Evershed effect causes the line cores to be blue-shifted on 
the limb-side of the penumbra and redshifted at the other 
side. The blue-shifted lines show a blue asymmetry that 
is negative, whereas the red-shifted lines show a positive 
asymmetry (Sanchez Almeida et al. 1996). While the con- 
tribution of the total light emitted is quite small, a spot 
does produce a hump in the profile of the photospheric 
line. Because we can measure the RV to a very high accu- 
racy, such effects can be detected. However, observations 
of the penumbra of spots close to the disk center show 
positive as well as negative asymmetry (Balasubramaniam 
1998). Thus it is not surprising that in the case of a star 
where there are numerous spots close to the pole, there is 
no correlation between the RV and the asymmetry. It is 
interesting to note that the bisector of the solar granula- 
tion shows the famous C-shape that will result in almost 
no shift between the upper half of the line and the lower 
half. We would call such lines symmetric. However, once 
the convective structures are spatially resolved, the red- 
shifted inter-granular lanes turn out to give rise to spectral 
lines with a negative asymmetry, and the granules give rise 
to lines with a slightly positive asymmetry (Guenther & 
Mattig 1991). 

Thus we presume that large and complicated flow pat- 
terns in the polar spots prevent us from observing a cor- 
relation between the RV and line asymmetry or depth of 
Ha. Nevertheless our observations show that the periodic 
RV signal can be detected if only the lines with positive 
asymmetry > 600 m/s arc taken into account. This indi- 
cates that there is a link between the asymmetry and RV, 
which implies that both are possibly caused by spots. 

Our observations also suggest that in the case of 
EK Dra it is impossible to account for the RV signal of the 
star spots to increase the sensitivity of a possible RV sig- 
nal of an orbiting planet because the correlation of either 
the depth of Ha or the line-asymmetry with the rotation 
phase is limited. 

6. Summary and conclusions 

We have studied the young star EK Dra spectroscopically 
and have resolved the binary system by means of speckle 
interferometry. This star is in fact a young star evolving 
analogously to the Sun, where the analogy in the evolution 
is even closer than previously thought. 

— By combining the speckle interferometry with old and 
new RV measurements we have for the first time a 
reliable mass estimate of EK Dra A and EK Dra B of 
0.9 ± 0.1 Mra and 0.5 ± 0.1 M^. 



— The orbit of the binary turns out to be highly eccentric 
with e = 0.82 ± 0.03. The period of the system is of 
45 ± 5 years, much longer than previously thought. 

— The high resolution spectra allow us to derive a new 
Teff of 5700 ± 70 K and a new value of log 5 of 
4.37 ± 0.10. Additionally, we measured a metallicity 
of [Fe/H] = -0.16 ± 0.07. The new values thus differ 
from those of Strassmeier & Rice (1998), who found 
Tcff = 5870 ± 50 K, log 5 = 4.5 and the metaUicity 
[Fe/H\ = 0.0. However, the spectroscopic distance de- 
rived using our new values now is in good agreement 
with the Hipparcos distance. 

— The Baraffe et al. (1998) models allow a mass deter- 
mination (1.O25±O.1OOM0) and an age determination 
of 35 Myr (30 Myr to MS) for the primary. 

— From our data we can exclude the presence of a third 
low-mass companion in the system as suggested by 
Metchev & Hillenbrand(2004). 

— We also see variations in the Ha line wings which lead 
to a temperature difference between the two spectra 
of 120 K. This temperature change could be caused by 
huge stellar spots. The effect on the line core and wings 
of the H/3 lines is not as strong. 

— The precise RV measurements taken in the course of 
the program shed more light on the activity of the star. 
We find that the primary star shows large amplitude 
RV variations which are presumably caused by its high 
activity level. In 2001 and 2002 the signal was periodic 
with a period of 2.767 ± 0.005 days. However, it seems 
that this period is slightly different from the photomet- 
ric rotation period (2.599 to 2.796, Strassmeier & Rice 
1998). This difference could be caused by differential 
rotation of the star. In 2002 this signal is significantly 
less obvious, and in 2003 almost absent. The periodic 
signal can be recovered from the whole dataset if only 
spectra in which the photospheric lines are asymmetric 
are used. We interpret this signal as being due to spots 
that are not located at the pole. However, the correla- 
tion between asymmetry and RV is weak. If the same 
analysis is performed on spectra in which the photo- 
spheric lines are nearly symmetric, we cannot recover 
the period signal of the RV variations. Also there is 
no obvious correlation between Ha and either the RV 
or the asynunctry. We interpret this by assuming that 
the star has large spots on the pole and possibly the 
Evershed-fiow of these spots affects the RV, as well as 
the asymmetry signal. We conclude that in the case of 
EK Dra the RV signal is related to the rotation pe- 
riod of the star and is caused by its activity cannot be 
simply subtracted from the RV signal to increase the 
sensibility for the detection of a possible third object. 
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